Conductance enhancement due to atomic potential fluctuations in graphene 
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We solve the Dirac equation, which describes charge massless chiral relativistic carriers in a two- 
dimensional graphene. We have identified and analysed a novel pseudospin-dependent scattering 
effect. We compute the tunneling conductance and generalize the analytical result in the presence of 
the tunable atomic potential of a graphene strip. The absence of back scattering in graphene is shown 
to be due to Berry's phase which corresponds to a sign change of the wave function under a spin 
rotation of a particle. We use the transfer matrix approach and find that the electric conductance 
of doped graphene increases due to atomic potential fluctuations. 
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INTRODUCTION 



Graphene is a single layer of carbon atoms densely 
packed in a honeycomb lattice. This material was found 
in its free state only recently, when individual graphene 
samples of a few microns in size were isolated by mi- 
cromechanical cleavage of graphite The current in- 
tense interest in graphene is driven by both the unusual 
physics involved and a realistic promise of device appli- 
cations. Among relativistic like phenomena observed in 
graphene so far, there are two new types of the integer 
quantum Hall effect and the presence of minimal metal- 
lic conductivity of about one conductivity quantum[2|. 
Transport properties of graphene are interesting because 
of their unique topological structure There have 

been some reports on experimental study of transport in 
graphene [7[. Various calculations have been performed 
to understand energy bands of graphene || Q. It has 
been successful in the study of various properties includ- 
ing the Aharonov-Bohm effect [10j, optical absorption 
spect ra II ill, quantum tunneling in graphene-based struc- 
tures fl2Hl4j , ferromagnetism [15| and instabilities in the 
presence and absence of a magnetic field [l6j . 

In a two-dimensional world, there are two basic mo- 
tions: forward and backward. Random scattering can 
cause them to mix, which leads to resistance. Just as 
we have learned from the basic traffic control, it would 
be much better if we could spatially filter the counter- 
flow directions. This effect could be used to produce a 
valley-polarized current out of an unpolarized stream of 
electrons due to line defect particularly in graphene []~Tl — 

The electronic transport properties of graphene with 
defects and impurities exhibit pronounced differences 
from those of conventional two dimensional electron sys- 
tems investigated in the past. The basic question that 
we seek to answer is what happens to the tunneling as 
we approach the Dirac point of zero carrier concentration 
with the presence of electrically charged scatters. In the 
present work we consider the effects of atomic potential 
Auctions in the setup depicted in fig. 1. At low doping 
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FIG. 1. Schematic of a graphene layer with a line tunable 
potential. The line potential is placed inside the graphene 
strip enhances the conductance in a vicinity of the Dirac point 
e <C hv/L, where L is the length of the strip with the width W. 
The indices A, B label the two sublattices of the honeycomb 
lattice of carbon atoms has different atomic potentials, which 
can be tuned. 



the conductance is determined by quasiparticle tunnel- 
ing, which is independent of the boundary conditions in 
the y-direction if L <C W. 



II. BACK SCATTERING IN GRAPHENE 

A. General consideration 

The quasiparticle excitation spectrum of the graphene- 
based junction consists of the positive eigenvalues of the 
Dirac equation. The Dirac equation has the form of two 
equations for electron u(y)a,b and hole wavefunctions 
v{t)a,b- In this study we restrict ourselves to the single- 
valley Dirac equation for graphene 



ihvF<yVu + Vu = eu 



(1) 



where u is a spinor of wave amplitudes for two nonequiv- 
alent sites of the honeycomb lattice. The Fermi energy e 
and the impurity potential V (x, y) in the graphene sam- 
ple (0 < x < L) are considered to be much smaller than 
the Fermi energy Ep in the ideal metallic leads (x < 



and x > L). For zero doping the conductance is deter- 
mined by the states at the Dirac point, e = 0. Transport 
properties at finite energies determine the conductance 
of doped graphene. 

Since for aspect ratios W L the boundary condi- 
tions in the y-direction are irrelevant, so we take peri- 
odic boundary conditions for simplicity. Different wave 
vectors g„=27rn/W (with n=0, ± 1, ± 2, . . . ) in the y- 
direction parallel to the interfaces(at the points x = and 
x = L) and not coupled, so we can consider each trans- 
verse mode separately. Without impurity in the bulk at 
a given energy e and transverse wave vector q we have 
up to two basis states 
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with the definitions 



a(e) = arcsin[fi,Wi?q/(e + Ep)], 
k(e) = (HvF)~ l (e + E F ) cos(a) 



(2) 



(3) 
(4) 



The angle a € (-7r/2,7r/2) is the angle between the ini- 
tial and final wave vectors = (k,q)\ x= o and k/ = 
(k, q)\x=L ■ With this sign convention the state u + move 
in the +x, while u~ move in the —x direction. The factor 
1 / y/cosa ensure that two states carry the same particle 
current. 

Considering a back scattering process k — > — k in an 
arbitrary external potential we will confine ourselves to 
states in the vicinity of the K point, but the extension 
to states near a K point is straightforward. We have 
a(k) = and a(— k) = n. Rotation in spin space can be 
obtained through the replacement (sj,kj) — ► (s/, — ky) 
corresponding to the electron motion of a time-reversal 
path. When the wave vector k is rotated in the anti- 
clockwise direction adiabatically as a function of time t 
around the origin for a time interval < t < T, the wave- 
function is changed into w s (k) exp(— i<J)), where u s (k) is 
the " spin" part of an eigenfunction of the Eq. Q] and <j) 
is Berry's phase. Choosing </>(k) = a(k)/2 in such a way 
that the wave function becomes continuous as a function 
of a(k)/2 and given by 



i I dt(u s (k(t)) 
'a 



d 
dt 



u.(k(t))) 



(5) 



This shows that the rotation in the k space by 2tt leads 
to the change in the phase by +n, i.e., a sign change. A 
similar phase change of back scattering process is the ori- 
gin of the so-called anti-localization effect in systems with 
strong spin-orbit scattering[27]. This anti-localization ef- 
fect was observed experimentally [28|, [2l| . The absence 
of back scattering can be destroyed by various effects. 
When the potential range becomes shorter than the lat- 
tice constant, the back scattering appears two reasons. 
The first is the appearance of inter- valley matrix elements 
between K and K' points [30l | . The second, which we in- 
vestigate in this work, is that the effective potential of an 
impurity for A and B sites in honeycomb lattice can be 



different. This causes mixing of the spin space and the 
momentum space. 



B. Long range impurities 

There is an absence of back scattering for long range 
impurities. This can be seen easily: using the formalism 
of the previous subsection we can simply show that the 
matrix element of the back-scattering process is zero 



-k|ff, mp |k) = -L 



dre 



Jk 



V(r) 
V(r) 



-»o k /2 
oiak/2 



c/2 



0(6) 



(here S is the unit cell area of the graphene lattice). 
This is rather dramatic since all real materials have im- 
purities, and it is normally these which determine the 
transport properties of the material. However, in the 
case of graphene the impurities don't scatter, i.e. they 
are effectively not there, and hence one has the possibility 
of coherent transport in graphene. 



III. TRANSFER MATRIX APPROACH 

We first determine the transfer matrix M n (x, 0) of the 
nth mode u n (x)e lqrlV through the undoped graphene rib- 
bon. The evolution of u n {x) inside the graphene sam- 
ple can be written as u n (x) = M n (x, 0)w„(0), where the 
transfer matrix M n (x) should satisfy a generalized uni- 
tarity condition 



M^ 1 = a x M^a x 



(7) 



In general Eq[6] (ignoring an irrelevant scalar phase fac- 
tor) restricts M to a three-parameter form 



(8) 



The real parameters a\, ct2, a depend on the boundary 
at the scale of the lattice constant and they cannot be 
determined from the Dirac equations. In this work we 
consider the ideal interfaces x = and x = L case, when 
a.\ = = oi2- From the Dirac EqQ]we obtain the differ- 
ential equation 



dx 



M n (x,0) 



is 
hv 



(9) 



with solution for V(x, y) = 

sin k n x ( ie 



M n (x, 0) = cos k n x ■ 



j^<Jx + q n o z ) (10) 



where k n — *J ((-^) 2 — ~q\ is the longitudinal wave vec- 
tor, q n is the transversal wave vector. The total trans- 
fer matrix through the impurity-free graphene ribbon 
(0 < x < L) is 



M tot = M(L,0)e iQ<T *M(0,0) 



(11) 
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FIG. 2. Gate voltage dependence of the conductance (calculated from Eq |17[) for the case of a smooth edge q n — {n + |) 
a.) for a different aspect ratios W/L — 5, W/L — 10, W/L — 20, and for incident angle a — 0°; b.), c), d.) for a fixed aspect 
ratio W/L — 10 and for a = 30°, a — 60°, a — 90° angles correspondingly. The horizontal asymptotes are indicated by dashed 
lines. 



Now we construct the transfer matrix for the electron 
(hole) quasiparticle excitations in the graphene ribbon for 
disorder potential case, where V is the charged impurity 
potential given by a diagonal matrix, i.e., 

V=( T*-*o) \ (12) 

where 7^4 (75) is the corresponding microscopic potential 
and smooth on atomic scales, which is localized along 
a line x — xq and placed in the zone corner K (K ) 
of sublattice A (£?). The Fermi-energy e and impurity 
potential V (Eq fT2|) in graphene sample (0 < x < L) are 
considered to be much smaller than the Fermi -energy Ep 
in the ideal metallic leads (x < and x > L). For zero 
doping the conductance is determined by the states at 
the Dirac point, e = 0. Transport properties at finite 
energies determine the conductance of doped graphene. 

The great physical and analytical advantage of the M 
matrix is that it is multiplicative. In order to obtain 
transfer matrix for disorder region (along line x — xq) 
graphene Md we integrate Eq. [1] 

rxo+u i-xq+v 

lim / (—ihvpa ■ yu + Vu) dx — lim / eudx 



Making use of well-known algebraic property for Pauli 
matrices 

°l = <4 =°l = -icr x crycr z = f J ° j (13) 

we easily obtain the final form for total transfer matrix 
M in disorder graphene sample 

M tot = M( J L,0)e^' Ta =M d e^ <Ta; M(0,0) (14) 

where 

Mi?) (l5) 

and ^a(b) — ~~i"tA(B)/ti v F- The transfer matrix of the 
whole sample is straightforwardly related to the matrices 
of transmission and reflection amplitudes 

M = M ™={%~l/f ) ( 16 ) 

The conductance of the graphene strip is expressed 
through the transmission amplitudes by the Landauer 



3 



formula 



IV. CONCLUSION 



N-l 

G = g J2 Tn > 9o=4e 2 /h (17) 

n=0 



where T = Tr[tH]. At the Dirac point V gate = the 
transmission probability for the case N 3> W/L reads 



(cosf +Psinf) 2 cosh 2 [7r(n+|)i] 



When parameters P = and a = the EqQ~8] reduces 
to equation (4), derived in the paper [31j . Parameter 
P =\ ja — Jb | reperesents atomic potential fluctuations, 
when finite P <C Hvf- 



In this paper we used the transfer matrix approach to 
calculate the tunneling conductance. We generalize the 
analytical result in the presence of the tunable atomic 
potential of a graphene strip and find that the electric 
conductance of doped graphene increases due to atomic 
potential fluctuations. The absence of back scattering in 
graphene is shown to be due to Berry's phase which corre- 
sponds to a sign change of the wave function under a spin 
rotation of a particle. We have identified and analysed 
a novel pseudospin-dependent scattering effect. The low 
cost experimental techniques allow for easy verification 
of the proposed hypothesis. 
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